Cross-Correlated Responses of Topological Superconductors and Superfluids 
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We study nontrivial responses of topological superconductors and superfluids to the tempera- 
ture gradient and rotation of the system. In two-dimensional gapped systems, the Streda formula 
for the electric Hall conductivity is generalized to the thermal Hall conductivity. Applying this 
formula to the Majorana surface states of three-dimensional topological superconductors predicts 
cross-correlated responses between the orbital angular momentum and thermal polarization (entropy 
polarization). These results can be naturally related to the gravitoelectromagnetism description of 
three-dimensional topological superconductors and superfluids, analogous to the topological mag- 
netoelectric effect in Z2 topological insulators. 

PACS numbers: 73.43.-f, 74.25.fc, 74.90.+n, 74.25.F- 



Introduction The quantum Hall effect (QHE) jl| is 
a prominent example of quantum phenomena charac- 
teristic of insulators with topologically nontrivial elec- 
tronic wave functions, a class of materials called topo- 
logical insulators (TIs) ||. In the QHE the Hall con- 
ductivity is quantized in units of e 2 /h at integer values 
equal to the topological number of bulk wave functions 
j| . The two-dimensional (2d) topological superconduc- 
tors (TSCs) and superfluids (TSFs) with chiral (p-wave) 
Cooper pairing arc superconductor analogues of the QHE 
and considered to be realized, e.g., in a thin film of 3 He 
A phase ||, [|, Sr 2 Ru04 ||, and the v — 5/2 fractional 
QHE fj. The topological nature of such TSCs and TSFs 
will manifest itself in thermal transport properties, such 
as quantization of the thermal Hall conductivity Q . 

Recent studies have shown that topological states ex- 
ist in time-reversal invariant and three-dimensional (3d) 
cases as well Q , and the systematic classification of them 
is established in terms of symmetries and dimensionality 
|J [n|. A key experimental signature of 3d- TIs is the 
topological magnctoelectric (ME) effect. Namely, the 
electromagnetic response of 3d- TIs is described by the 
axion electrodynamics ^ [0], p^[ , 
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with 6 = ir, the possible nonzero value in time-reversal in- 
variant systems (mod 2ir). The effective action ([!]) leads 
to the surface QHE that induces the topological ME ef- 
fect as M = (e 2 /2hc)E and P = (e 2 /2hc)B where M 
and P are the magnetization and electric polarization, 
respectively. 

An example of 3d-TSFs |l3| is the B phase of super- 
fluid 3 Hc B. In addition, the newly found superconduct- 



ing phase in Cu-doped Bi2Se3 |14 has been proposed to 
be a 3d-TSC @. The topological nature of 3d-TSCs 
manifests itself in their coupling with the gravitational 



field, which is described by the term similar to Eq. (Q), 
the gravitational instanton term [[l6| [l7] . 

In this work we shall develop the response theory of 
TSCs, which reveals the cross-correlation between ther- 
mal and mechanical (rotational) responses, and is the 
direct equivalent of the topological ME effect in TIs. We 
first generalize the Streda formula |Q to the thermal 
Hall conductivity in 2d systems. By applying it to the 
surface states of a 3d-TSC, the cross-correlated responses 
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FIG. 1. (color online) Electronic responses in (a) two- 
dimensional (2d) and (c) three-dimensional (3d) topological 
insulators (TIs) (left) and thermal and mechanical (rotat- 
ing) responses in (b) 2d and (d) 3d topological superconduc- 
tors (TSCs) (right). In (b), the arrow along the boundary 
indicates the chiral Majorana edge channel with dispersion 
E = vpe- In the 3d-TSC (d), temperature gradient induces 
surface thermal Hall current je- A uniform mass gap is in- 
duced in the surface fermion spectrum by doping magnetic 
impurities near the surface of the 3d-TI and TSC such that 
spins are all pointing out or in (red arrows). 
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TABLE I. Comparison between cross-correlation in topolog- 
ical insulators (TIs) and topological superconductors (TSCs) 
in two (2d) and three spatial dimensions (3d). In TSCs the 
orbital angular momentum L and entropy S (thermal po- 
larization Pe in 3d) are generated by temperature gradient 



(Eg 



-T 1 VT) and by rotating the system with angular 



velocity fi a . In analogy with the magnetoelectric polarizabil- 
ity Xe b m 3d TI gravitomagnetoelectric polarizability Xe b g can 
be introduced in the 3d TSC (Right-bottom). Note that the 
relations for TSCs hold also for the thermal response of TIs. 
These responses are characterized by the topological integers 
and quantized. See Eq. (g) for 2d-TSCs and Eq. @ for the 
surface of 3d-TSCs 



of 3d-TSCs arc identified with those between the orbital 
angular momentum and the thermal (entropy) polariza- 
tion, generated by temperature gradient VT and by ro- 
tating the system with angular velocity f2, respectively, 
as shown in Fig. [l]. Our main findings are summarized 
in Table |IJ Their derivations will be given below. 

Thermal Hall conductivity Relation between the Hall 
conductivity a H and the magnetization M in the quan- 
tum Hall regime is known as the Streda formula |lq] : 



dM z 



(2) 



where fi is the chemical potential, e (< 0) is the electric 
charge, and c is the speed of light. The magnetization 
is evaluated by M = (e/2c)Tr[£>(/i — TL)x x v], where 
v = (i/h)[H,x], TL is the single-particle Hamiltonian, 
and 9(x) is the step function. The relation (0) can be 
understood by identifying the Hall current j = a H E x z 
with j = cV x M = -cdM/d^i x V/i. 

For the thermal conductivity k h = j'^/d y T, we will 
show below that a similar relation holds: 



dM z 
2 dT 



(3) 



The thermal current is j T = j E — (n/e)j, where j is 
the electric charge current. The energy-current j E is 
defined as j E = cT a0 , where T a0 is a spatio-temporal 
component of the energy-momentum tensor, T^ v = 
(8 C/ 58^)3"^ — g^ v £, satisfying the continuity equa- 
tion, d^T^ v = 0; £[^,8^] is the Lagrangian density 
of the system. The moment of the thermal current in 
Eq. (pD is defined as M T = M E - (fj,/e)M, where 



to the charge Hall conductivity (||) , the average has to be 
taken at finite temperature: (• • • ) = y, n f(e n ){n\ ■ ■ • |n) 
where e n and \n) are eigenvalue and eigenstate of the 
Hamiltonian TL, and /(e„) is the Fermi distribution func- 
tion. Below, we will consider the part of the thermal cur- 
rent carried by je, specializing to the case of [i = 0. For 
insulators, this means we count the total number N of 
particles such that when the chemical potential is within 
a gap, N — 0. In the Bogoliubov-de Genne theory of 
superconductors, due to particle-hole symmetry, this is 
always true, and hence jr = 3e and Mt = Me- 

Let us now introduce a gravitomagnetic field B g pj| 
which is conjugate to Me so that the variation of free 
energy is dF = -SdT - M E ■ dB n . Equation (0) is 
written as k h = c(dM E /dT) B * 



dB g . Equation ( 
= c(dS/dB* q ) T , or 
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by introducing dQ = TdS and d<p = dT/T, where <fi is 
a fictitious gravitational potential that couples to ther- 
mal energy density Q [p0| . Equation (|J) is the thermal 
analogue of the Streda formula for the charge Hall con- 
ductivity, a H = ecdM z /dfi = ecdN/dB z , in that Q is 
the zeroth component of the energy current as eN is in 
the charge current. 

To see the physical meaning of B g and Me, note that 
the definition of is similar to the orbital angular mo- 
mentum: = (1/c) (x^T 0v - x u T°' J ') . Indeed, when 
there is a rclativistic invariance, the energy-momentum 
tensor can be symmetrized so that = T vfl , and 
thus M*" = (c/2)L^ u . While this is not the case in 
condensed-matter systems in general, in so-called pseudo- 
rclativistic systems where the Lorentz invariance is real- 
ized at low energies, electrons or quasiparticles obey the 
Dirac or Majorana equation. In these systems the Fermi 
velocity v plays a role of the speed of light c, and the 
tensor is related to the orbital angular momentum 
as Me = (v/2)L, and B g can be understood as the an- 
gular velocity vector of rotating systems, B g = [2/v)H. 
For a system rotating with the frequency f2 = fl z z, this 
can also be understood by making a coordinate transfor- 
mation from the rest frame to the rotating frame in which 
the metric in the polar coordinates (t, r, ip) takes the form 
ds 2 ~ v 2 dt 2 — 2Q z r 2 d<pdt — r 2 dtp 2 — dr 2 . One can then 
read off, from the definition of the gravitoelectromagnetic 
field, the non-zero gravito gauge potential = Q z r/v 
Jl9| . In Cartesian coordinates, A g = (l/v)fl z z x x, and 
B g = 
written as 



V x A g = (2/v)VL z i. Consequently, kh can be 



i ^ - iT ) , and MJ 



M\ 2 . 



In contrast 



This is one of the main results of this work. 

Although it is necessary to have (pseudo) Lorentz in- 
variance to identify, at the operator level, Me with the 
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angular momentum L, the relation (Q) can in fact be 
derived for a wider range of systems, e.g., by assuming 
that the edge state in the disk geometry is described by 
chiral conformal field theory This indicates the va- 
lidity of Eq. (§) in many 2d topological phases, a rep- 
resentative example of which is the 2d chiral p-wave SC 
with l z = — 1 pairing as shown in Fig. 0(b). Near the 
edge, there exist chiral Majorana modes M described by 
-Hedge = (1/2) Jq dxip(— ihvd x )ip, where ip is a single- 
component real fermionic field and L is the circumference 
of the edge. Since T m = T w = T 01 = (-ifw/2)ipd x ip, 



( JE ) = (v/L)(H, 
leading to 



) = (v 2 /2)(L^ e ) 



dT 



^ 2 k 2 B T 2 /12h, 
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6h 
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which is a half of the quantized value for conventional 
chiral fermions. Here, £^ dgG is the edge modes contri- 
bution to the total orbital angular momentum per unit 
area. 

Majorana fermion on the surface of 3d-TSC/TSF 
Let us now derive Eqs. (^jS]) from a microscopic the- 
ory via Kubo formula. We work with an example, 
the 2d massive Majorana fermion system realized on 
the surface of 3d-TSCs fi ||, |22|], anticipating to ap- 
ply (@-||) to derive the cross-correlated responses of the 
3d-TSCs. It is described by H = (1/2) / d 2 xiP T Hip, 
where ip T = (ijjf, 4>i) is the real spinor field satisfying 
{ij) a (x),il>p(x')} = 6 a/3 S(x - x'), and 



H = —ihv(o z d x + cr x dy) + ma v 



(7) 



The mass to is due to the interaction with magnetic fields 
or magnetic moments perpendicular to the surface [Fig. 
|l|(d)] and thus breaks time-reversal symmetry. To study 
thermal transport, we introduce coupling with the ficti- 
tious gravitational potential (j) in the Lagrangian ^3| : 
C = (l/2)ip T [ihd t - U- (1/2){0,^}]V>. The energy- 
current operator is given by (a, b = x,y) 

+ ip T - [H{v a x b + 3x b v a ) + H.c.l tpdhcf). (8) 
8 

The first term is non- vanishing even in the absence of 
gravitational potential; the second and the third terms 
are proportional to V</>, which are analogous to the dia- 
magnctic charge current in the presence of a magnetic 
field. To evaluate kh = —(j%)/(Td y (j)), we apply the 
Kubo linear response formula to the first term, while the 
second and third terms need to be averaged in thermal 
equilibrium plj . When the mass gap is large enough, the 
density of states of Majorana fermions vanishes at e = 0, 
and the thermal Hall conductivity is given by kh = 
vdM^/dT, where M% = (l/4«) £„ f(s n )s n (n\{xvv - 
yv x )\n), and (x\n) = u n (x) is the exact eigenstate of 
the Majorana Hamiltonian (^), Hu n (x) = e n u n (x). 



Quite generally, one can derive, in the limit T — > 0, 
the relation 



K H 



hir 2 k%T sr ^ 2lm[(n\v x \m)(m\v v \n)} 



6L 2 



(9) 



where L 2 is the of the surface. Apart from 

the factor 7r 2 fc^T/6, the right hand side resembles the 
Kubo formula for the electrical Hall conductivity, which, 
however, is not a well-defined quantity for Majorana 
fermions; nevertheless Eq. @ can be regarded as the 
generalized Wiedemann-Franz law to Majorana fermions 
j| ||, |2^, p3[ . Compared to the electron systems, there 
is an extra factor of 1/2 due to Majorana nature. Since 
oh of the massive Dirac fermion is a H = sgn(m)e /(2h), 
Eq. (^) immediately gives 



, H = sgn(m) y -T 



(10) 



for the massive Majorana fermion. There is a factor 1/2 
compared to the 2d result Eq. (||). 

Cross-correlated response of 3d-TSC/TSF Let us il- 
lustrate physical implications of Eq. (^) by studying the 
responses of 3d-TSCs to the temperature gradient and 
the rotation. For simplicity, we consider a sample with 
cylindrical geometry with height £ and radius r as de- 
picted in Fig. 0(d). We assume that magnetic impuri- 
ties are doped near the surface and their spins are all 
pointing out or in so that uniform mass gap is formed on 
the surface. Let us first introduce the temperature gradi- 
ent in the z-direction, which generates the energy-current 
3e = Kfid z T on the surface. Since Je/v 2 corresponds to 
the momentum per unit area, total momentum due to the 
surface energy-current is P v = (2irrl)jE /v 2 and thus the 
induced orbital angular momentum per volume is given 
by 



rP w 



S2 = 



d z T. 



(11) 



Similarly, upon rotating the cylinder with f2 = ft z z 
(without temperature gradient), applying Eqs. (||) and 
(||) to the top and bottom surfaces, we obtain the induced 
thermal energy density (the induced entropy change) lo- 
calized on the top and bottom surfaces, 



AQ(z)\ 



where k\j 



2Tfl z 

[n^Siz - 1/2) + K h H S{z + £/2)} , 

v 

(12) 

-Kj^ as the spins on the top and bottom 
surfaces are pointing to the opposite directions (different 
signs of to); sec Fig. |l|(d) |24| . 

In terms of the gravitoelectric field E g = — T^ 1 VT and 
the momentum of the energy current Me, Eq. ( |Tl| ) can 
be written as Me = {Tkh /v)E g . Further introducing 
the thermal polarization Pe by AQ = -V-P^, Eq. (TIT 



can be written similarly as Pe — {Tkh /v)B g . These 
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highlight the correspondence between TIs and TSCs, 
TI: = ^ * TSC: ^ = 9. (13) 



dE b dB b 



dE b dB b 



Since the orbital angular momentum is given from the 
internal energy functional by L a = ~8Ug / 8Vt a , the cou- 
pling energy of the temperature gradient and rotation 
velocity is written as 



k 2 T 2 
-^—6E„ 

24hv 9 



This is analogous to Eq. (|l|) with e 2 /he O (nk B T) 2 /6hv 
and 9 = ±n playing the same role as in Eq. (0). In 2d- 
TSC cases, the corresponding term is written as t/xsc = 
J d 2 x(2/v 2 )TK H (f>H z . This is the thermodynamical ana- 
logue of the Chern- Simons term 0, pi B5) . 

Angular momentum paradox in 2d-TSC/TSF These 
terms can be related to the problem of total angu- 
lar momentum of the ground state of chiral supcrflu- 
ids. For simplicity we consider the 2d case. The an- 
gular momentum per unit area of the ground state of 
the chiral p-wave state has been proposed to behave as 
L Z (T = 0) = -(hn/2)(k B T c /E F )~< with different expo- 
nents 7 = 2, 1 and U, where n is the number of par- 
ticles. The controversy of 7 has not yet resolved. By 
integrating Eq. (||) from to T, and using A = HvkF, 
one obtains L Z (T) -L z (0) = (nhk 2 F /6)(k B T/ A) 2 , at low 
temperature, which is consistent with a numerical simu- 
lation |26| . If we extrapolate this relation to the critical 
temperature T c at which L Z (T C ) vanishes, we obtain the 
angular momentum of the ground state, up to a numeri- 
cal factor, L Z (T = 0) h(irk 2 F /2)(k B T c / A) 2 hn/2, 

indicative of 7 = 0. (Here we should keep in mind that 
we used kh oc T which is, however, valid only at low 
temperature T <C T c .) Since the chiral Majorana edge 
modes of a 2d-TSC are moving in the opposite direc- 
tion to the rotation of Cooper pairs [Fig. 0(b)], at finite 
temperature, they contribute to reduce the total angular 
momentum from its ground state value at T = 0. 

Possible experiments We conclude by discussing pos- 
sible experiments to probe the cross-correlated response 
of TSCs. For the 2d case, let us assume a 2d-TSC is 
rotating with frequency Vl z . As il z increases by Ail 2 , 
we predict the temperature changes as AT = AQ/C = 
{2khT / Cv 2 )AVL z 1 where C is the heat capacity of the 
2d-TSC sample. As the heat capacity of a fully gapped 
superconductor is small |j) , this temperature change may 
not be so difficult to measure. 

For 3d-TSCs, the Einstein-de Hass effect will reveal a 
relationship between magnetism and angular momentum 
of the system, as it has been used for the study of fer- 
romagnetic materials. We assume a cylindrical 3d-TSC 
suspended by a thin string and apply thermal gradient 



[Fig. 0(d)]. This induces surface energy current with an- 
gular momentum L z , according to Eq. (|ll]) . By the con- 
servation law of total angular momentum, it must be 
compensated by a mechanical angular momentum of the 
material, which can be directly measured in principle. 
Its inverse effect is the generation of thermal polariza- 
tion (entropy polarization) by rotating the system. To 
observe the effect, the frequency il z should be lower than 
the critical frequency f2 c i above which vortices are intro- 
duced in the bulk of the sample that will generate the 
energy current in the z direction. The key point is that 
magnetization has to be induced properly on the surface, 



(14) all pointing out or in 24 



In summary, we have found nontrivial correlated re- 
sponses of TSCs to thermal gradient and mechanical ro- 
tation which arc analogous to the topological ME ef- 
fect. We have proposed possible experiments which probe 
these topological responses. 

While we were finalizing the paper, a preprint |27|| 
appeared in which the energy magnetization Me is 
discussed while its relation to the angular momentum 
L and TSCs arc not discussed. This work is sup- 
ported by MEXT Grand-in-Aid No. 20740167, 19048008, 
19048015, 21244053, Strategic International Cooperative 
Program (Joint Research Type) from Japan Science and 
Technology Agency, and by the Japan Society for the 
Promotion of Science (JSPS) through its "Funding Pro- 
gram for World-Leading Innovative R&D on Science and 
Technology (FIRST Program)." 
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THERMAL STREDA FORMULA FROM THE 
CFT PARTITION FUNCTION 

In this appendix, we will derive the thermal Streda 
formula (||) of 2d topological liquid by making use of 
the partition function of the chiral conformal field theory 
(CFT) on the edge. Our derivation closely follows the 
calculations of the thermal conductivity, the Coulomb 
blockade current peak, and the thermopower by Cappclli 
et al 

While the Streda formula can be derived solely from 
the Kubo formula in the bulk, as illustrated in the text 
in terms of the 2d Majorana fermion system, the purpose 
of such complementary derivation is two-fold: (i) In the 
bulk of a topological fluid, the dispersion of gapped quasi- 
particle excitations is not universal and one can adiabat- 
ically deform it without changing the topological proper- 
ties of the system. On the other hand, the dispersion of 
gaplcss excitations on the edge is, at low energies, quite 
universal with emergent Lorcntz invariance. 

(ii) While the derivation in terms of the Kubo formula 
mainly focuses, at least in the form presented in the main 
text, the situation where there is no interaction among 
quasiparticles, the derivation in term of the edge theory 
applies to an arbitrary (chiral) topological liquid, includ- 
ing strongly interacting ones such as the fractional quan- 
tum Hall liquid. 

For the Streda formula for the electrical Hall conduc- 
tivity in the quantum Hall systems, and for the spin Hall 
conductivity in the quantum spin Hall systems, an ap- 
proach similar to the calculations we present in this sec- 
tion can be developed. Below, we will set H = c = fcs = 
1. 



partition function 

Let us consider a topological fluid in the annulus geom- 
etry which is boundary by two (inner and outer) edges. 
The low-energy physics is governed by the two edge 
modes. The Hamiltonian for the edge modes may be 
written as 



H 



vi t c 



V2 
R-2 



24 



(15) 



where Lq and Lq represents the energy-momentum ten- 
sor for the outer and inner edges, respectively; v\ and vi 
are the fermi velocity for the outer and inner edges, re- 
spectively; R\ and i?2 are the circumference of the outer 
and inner edges, respectively; By "adjusting" the veloc- 
ities of propagation of excitations, we may set v\/ R\ = 
V2/R2 = v/R. 

We now consider the grand canonical partition func- 
tion for the edges defined by 



z _ rj, re -i3[H-n z L z -ti{N+N)-V (N-N)] 



(16) 



where we have included the frequency f2 z , the chemical 
potential [i, and the electro potential difference between 
edges V ; they couple to the total orbital angular momen- 
tum, the total electron number, and the difference be- 
tween the electron numbers in the inner and outer edges; 
N and N represent the total electric charge of the inner 
and out edges, respectively. 

Below, we will organize these intensive variables in 
terms of two complex numbers r and 



v + in z /(irR) 
R ; 



2mC=-P(V a +ifi), (17) 



where we have "analytically continued" il z — > i£l z , and 
fx — > i/j.. The grand partition function Z(t, £) is then 
given by 



Z(r,C)=tr 



3 2iriT(L -c/24)-27rir(L -c/24)+CA r +CA r 



(18) 



In general, for a rational topological phase, where the 
number of types of bulk quasiparticles is finite, the an- 
nulus partition function can be computed as 



(19) 



where XaC^C) i s a character associated with the quasi- 
particle of type a, 777 is complex conjugation, ■ ■ - c rep- 
resents particle-hole conjugation, n a g is a non-negative 
integer, and the summation runs over a finite set of quasi- 
particle types. More specifically, the character x (t, Q is 
given by 



D 2mT(L -c/24)+CN 



(20) 
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where tr a [■ • • ] represents the trace within the Hilbcrt 
space associated to the particle type a. 

On the disk geometry bounded by only one edge, the 
partition function does not consists of holomorphic and 
antiholomorphic parts, but reduces to 



Z(T,0=Xa(r,0, 



(21) 



where the label a depends on quasiparticles present in 
the bulk. In particular, in the absence of quasiparticles, 
we have a = which represents the identity or "vacuum" 
particle. 



Streda formula 

In order to compute the thermal Hall conductivity, 
Kh = BJe/QT, we focus on the case with single edge 
(i.e., disk geometry), and compute the expectation value 
of the energy current Je, which can be written as 



Je = ve 



(22) 



where e is the average chiral energy density on a single 
edge and v is the velocity of the edge excitations. This 
can be computed from the chiral partition function ( |2l| ) 

as 



d 



In Xa, 



2irRfi dr] 
where we have introduced rj 

r] := ifl z /(TrR) 



(23) 



(24) 



for notational convenience. (Henceforth we set ( = ( 
0). The thermal Hall conductivity is then given by 



d d 1 



o ln Xa- 



2ttR dT drj /3 
On the other hand, the entropy is give by 

S = — In Z = In y„ . 

dT /3 dT (3 A 



(25) 



(26) 



Taking the derivative with respect to il z on the both 
sides, 



v 2- JUL 

2 dn z 



v 2 i d d 1 



lnXa- 



(27) 



We thus obtain the Streda formula for the thermal con- 
ductivity: 



d,J E 



dT 



f_ dS_ 
Y dip 



(28) 



We now evaluate the thermal conductivity in the large 
size R — > oo and in the low-temperature /? — > oo limit. 



It is important to take the R — > oo limit first, so we take 
the limit 



v/3/R < 1 



(29) 



first and then take the small T limit. The limit v/3/R — > 
can be evaluated by making use of the modular transfor- 
mation of the characters 



X 



(30) 



where is the modular S- matrix. The vf3/R — > limit 
of the RHS can be evaluated as 



Xa(q) = tr 



Lo-c/24 



ln 9 + ln(a;0|g Lo |a;0) 



24 
( _ ^4 + 



(31) 



where we assume the modulus of the complex parameter 
is small, |g| —> 0, and |a;0) is the highest 



= p -(2 7 r l /r) 



q := e 

weight state, and Lo| a ;0) = h a \a;0). Thus, 
Xa(r)^J2 S " e ~ {2WiM{ ~^ +hb) ' 



(32) 



In the limit Rc(27ri/r) = (4?r 2 R) / (v /3) -> oo, b = (the 
identity particle) dominates in the summation (unless 
S° a ^0). Thus, 



Extracting the temperature dependence, 



c d 2iri 2irT 2 



e ~ — - 



2tt.R24<97 t 
Finally, kh is given by 



v 24 



d 



d 2ttT 2 c 



ttT 



dT v 24 6 
Alternatively, one can obtain the entropy first: 

_^ 9 \ c 2iri ^ c 
~ T lfr) 2i~ ~ 24~' 



(33) 



(34) 



(35) 



(36) 



where, again, we take the limit v/3/R — > 0. The entropy 
can be expanded in small rj as 



6 p V 



We thus conclude 



u^_l dS 
in agreement with Eq. (|35|). 



C7T 

H 2 ttR dirj 6 



(37) 



(38) 
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THERMAL STREDA FORMULA FROM THE 2D 
MAJONARA HAMILTONIAN 

Streda formula 

The energy-current operator is divided into the follow- 
ing two terms: (a, b = x,y) 



where 



Je (o) 

-a 

Je(i) 



Je — Je(o) + Jjs(i)) 



*l> T -{v a ,H}i> 



(39) 



(40) 



il> T - [H{v a x b + 3x b v a ) + H.c.] i/>db(f>,(41) 



and ?j a = (i/h)[H, x a ]. ^ is non- vanishing even in the 
absence of gravitational potential, while ^ is propor- 
tional to V(f>. To evaluate kh = — {Je) I {T d y 4>) , we apply 
the Kubo linear response formula to j% ^ : 



(o) ) 
d b (j) 



i v f(E n ) - f(E m ) (E n + E„ 

nm 

(n\v a \m)(m\v b \n) 



E n - E m + ir) 



(42) 



Here (x\n) = u n (x) is the exact eigenstate of the Ma- 
jorana Hamiltonian, Hu n (x) — e n u n (x). On the other 
hand, j^^i) needs to be averaged in thermal equilibrium: 



4^^ /(£; " ) T (n|[wa ' ub]l " > ^- (43) 



One obtains the thermal Hall conductivity 



and G±(C) = (C — T~L ± irj)" 1 . It is convenient to write 
the thermal conductivity as a sum of the following two 
terms: 





1 


K H — 


2T 




d 


K H — 


df 



df(e) e 2 
de ~2 
1. 



B(e), 



cfe/(e)-Tr e(xu y - yv*)^ - H) 



where 



dM% 



(46) 



(47) 



M £ = ^£/( e »)e»H0«' 1 ' -2/^)1")- (48) 



At low temperature, the Fermi distribution function can 
be expanded as follows: 

f(Sn - C) = 0(C - En) - y T2 ^(C - En) (49) 



and thus 



^ 2 fc 2 T 2 



G 



-fl(C = o). 



(50) 



(j E (i))o = ^^f(E n )E n (n\(x a v b - x b v a )\n)d b <t> = f 



When the mass gap is large enough, the density of states 
of Majorana fermions and Eq. ( [45[ ) vanish at e = 0, and 
hence = k|J while = 0. 



Generalized Wiedemann-Franz law 

Here we derive Eq. (^|). Noting the identity 

v y S(C - H) - H.c. 



— OO 



dC ^Tr 



(n\v x \m) (m\v y \n) 



(e 



zr/) 2 



- C. 



(51) 



in Eq. and using Eq. (H|), one obtains Eq. (§). 



de 



-3/(e)\ e 



.2 f£ 



where g | f§ 



IdB 1 



(xz/ y — yv x ) 



d( A(C), (44) 



g(C - H) 



5(C) = -^Tr{KG+(CK " «*G-(C)«1*(C - W)}, 



(45) 
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